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The suppression of single inclusive hadron spectra in heavy-ion collision as compared to the 
scaled expectation from proton-proton collisions has long been regarded as an interesting tool to 
study properties of the bulk matter in heavy-ion collisions. However, the limitations of this class 
of observables has become increasingly obvious, and both experimental and theoretical efforts are 
now made to go beyond single hadrons to fully reconstructedjets. Monte-Carlo (MC) simulations 
of in-medium parton showers are currently considered the most promising tool to theoretically 
access jet physics in heavy-ion collisions. In this paper, I review some of the first results obtained 
with the MC code YaJEM (Yet another Jet Energy-loss Model) for both single hadron and jet 
observables. 
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1. Introduction 

Since the suppression of high transverse momentum (Pj) hadrons in heavy-ion (A-A) colli- 
sions as compared to the scaled expectation from proton-proton (p-p) collisions has been suggested 
theoretically as a tool to study properties of the hot and dense QCD matter produced in A-A colli- 
sions [jj]], a substantial effort has been made on the experimental side to measure high Pj hadrons. 
Current A-A data include the suppression of single inclusive hadrons as a function of Pj, both 
averaged [Q] and as a function of the angle of the hard hadron with the reaction plane and the 
suppression of back-to-back angular correlations 

Single hadron observables and back-to-back correlations are well described in detailed model 
calculations using the concept of energy loss ^ §,01, i.e. under the assumption that the process can 
be described by a medium-induced shift of the leading parton energy by an amount AE, followed 
by a fragmentation process with the vacuum fragmentation function of a parton with the reduced 
energy. However, there are also calculations for these observables in which the evolution of the 



whole in-medium parton shower is followed in an analytic way ffl, \A 1C]. A recent comparison 



study [11] demonstrates both the success of different models in describing the available data on 
single hadron suppression and the limitation of these observables in terms of constraining properties 
of the medium, as models with quite different physics assumptions are able to account well for the 
available data. 

One way of resolving this situation is to move from measurements of single hadron spectra 
to fully reconstructed jets. While the former measurement is sensitive primarily to the leading 
hadron of the jet from which energy is lost, in the latter measurement also subleading hadrons 
are observed and hence the redistribution of energy by the medium from the leading hadron to 
subleading hadrons or into the medium can be traced. On the experimental side, first measurements 
of jets have become available [|l2|]. MC simulations of in-medium parton showers are considered 
the appropriate tool to discuss this physics for a number of reasons. One important argument 
is that explicit conservation of energy and momentum may be a more important factor for jet 
observables than quantum coherence. Energy momentum conservation is not usually present in 
analytic calculations as they are commonly done in an asymptotic limit, but is easily accounted 
for in a MC code. On the other hand, quantum coherence is more easily accounted for in an 
analytical calculation than in a MC code [|13|]. Another crucial factor is that the experimental 
procedures to reconstruct jets in an heavy-ion environment can be readily applied to the output to 
a MC simulation, but it is currently unclear how to compare data with analyical calculations in 
an unbiased way. Thus, at present a number of MC codes to simulate jets in medium are being 
developed, among them JEWEL [§J]], YaJEM [0, f5|, and Q-PYTHIA [0. In this paper, I 
discuss the physics underlying YaJEM and the modifications to jets in heavy-ion collisions based 
on the results of the code while making the connection to older results based on the energy loss 
picture. 

2. Parton shower and fragmentation in vacuum 

In factorized Quantum Chromodynamics (QCD), the production of hadrons from a p-p colli- 
sion can in leading order be written schematically as 
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do NN^h+x = ^f^ixuQ^^fj/^^^aij^f+k^Df^Qj). (2.1) 

m 

Here, fi(j)/N( x i(2),Q 2 ) stand for the initial distributions of partem species i,j with momentum 
faction jqp) in the nucleon probed at scale Q 2 . Cij^f + k is the perturbatively calculable hard process 
of a reaction in which partons i,j scatter into partons f,k and D v "^ h (z,Q 2 ) is the distribution of 
hadrons h with momentum fraction z produced from a parton / probed at a scale \i 2 . Both / and D 
encode non-perturbative input, in particular the physics of hadronization and confinement, but their 
scale evolution with Q 2 or Q 2 is perturbatively calculable. Or, in other words, encodes both a 
perturbative partonic evolution from a parton / at an initial high virtuality scale Q 2 into a shower of 
partons at a lower scale Qq where perturbative QCD ceases to be applicable and a non-perturbative 
evolution of a parton shower at scale Qq into a hadron shower where the relevant scale is finally set 
by the squared hadron mass mf. 

If one wants to study the effect of the hot QCD medium produced in A-A collisions on this 
expression, it is useful to analyze these virtuality scales. By the uncertainty relation, the relevant 
timescale at which a fluctuation at a scale Q develops is l/Q in its own restframe, in the lab 
frame an additional Lorentz factor E/Q appears where E is the energy of the parton (hadron), thus 
for the evolution time X ~ E/Q 2 . Inserting typical energies E of 10-20 GeV for hard probes at 
RHIC and an initial virtuality scale <2, of the same order which is relevant for the hard process 
itself, one finds timescales T< 1 fm/c, i.e. the hard process occurs long before the thermalization 
timescale of x ~ 0.6 — 1 fm/c at which the medium is produced. On the other hand, inserting the 
pion mass as a typical hadronic scale into the expression, Q~m n , one finds x 3> x me d 10 — 20 
fm/c, i.e. the hadronization timescale is much greater than the lifetime of the medium. Thus, it 
follows that neither the parton distributions (which are an initial state effect relevant before the hard 
process), nor the hard process itself nor hadronization are affected by the medium, but chiefly the 
partonic evolution from an initial scale <2, down to the hadronization scale Qq. It is this evolution a 
computation of medium effects needs to focus on. 

Parton showers in vacuum can be modelled as a series of branchings a —> b,c with the three 
reactions q — > qg,g — > gg and g — > qq. Conveniently, the evolution is discussed in terms of two 
evolution variables in momentum space, t = InQ 2 / 'Aqcd as a measure of the (decreasing) virtuality 
scale and z the energy (momentum) splitting variable in each branching as E\, = zE a and E c = 
(1 —z)E a . In these variables, the differential branching probability at scale t is then given by 

dPa = £ ^P a ^bc(z)dtdz (2.2) 
b,c 171 

with the splitting kernels calculable in perturbative QCD as 

Pq^qgiz) = - -j— J Pg-+gg(z) = 3 j Pg^qqiz) = — {? + (1 - z) ). (2.3) 

The probability density that, given an initial scale tj n , the next branching occurs at a lower 
scale t m is then given by 



dPa 

dt m 





exp 


b,c 





b,c 



(2.4) 
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where / a ->6c (tm) is Eq. ( |2l2| ) integrated over z subject to kinematical bounds given by the parton 
virtualities. This expression has two terms, the first one stands for the actuai branching probability 
at t m whereas the second one, the so-caiied Sudakov formfactor, stands for the probability that 
no branching has occurred between ?,„ and t m . Eqs. (2.2,2.3,2.4) are solved in order to find the 
evolution of a parton shower in vacuum, and they need to be modified in order to compute the 
effects of the medium on the shower. If the expressions are solved in a MC simulation, quantum 
coherence can at least be approximately be included by imposing an angular ordering condition 



Zb{l-Zb) 1 -Za 



Ml 



> 



ZaMl 



(2.5) 



at each branching. YaJEM uses the PYSHOW routine Q180 which is part of the PYTHIA 



package [ |19| ] to solve the evolution equations in vacuum and a modifications of this routine to 
solve them in the medium. 



3. Energy loss as an approximation to in-medium shower evolution 

If one is interested in computing the medium-modification of single inclusive hadron spectra 
or back-to-back correlations, one can make an approximation in computing the medium-modified 
fragmentation function. The reason is that the primary parton momentum spectrum is steeply 
falling with parton momentum pj. Therefore, hadron spectra are dominated by comparatively 
hard fragmentation, for fragmenting quarks at RHIC kinematics z ~ 0.7. This means that one is 
dominated by events in which a single parton carries most of the momentum of the shower, and 
the fragmentation function chiefly parametrizes the hadronization of this leading parton and not 
so much the partonic shower. Any medium effect can then be pictured as energy loss from this 
leading parton. This is in contrast to an unbiased hard event in which the shower momentum 
is shared among many partons, thus leading to multiple low Pt hadron production. In the latter 
situation, clearly the assumption of energy loss from a leading parton is inadequate. 

In the energy loss approximation, the key quantity characterizing the medium is the geometry- 
averaged energy loss probability density (Pf(AE,E))T AA which appears in the factorized computa- 
tion between hard process and fragmentation. The medium-modified production of hadrons can 
then be computed (schematically) from 

d<r h+X = Zda^r f+X ® (Pf(AE,E)) TM ® D}Z h (z, /i 2 F ) (3.1) 
/ 

and 

do^ f+X = E/i/A(*i , Q 2 ) ®fj/ A (x2,Q 2 ) ® o iHf+k . (3.2) 

ijk 

The nuclear suppression factor Raa can in this framework be written as 

RAAiPr) = da^ h+x (P T )/da NN ^ h+x (P T ). (3.3) 

The averaging of the energy loss probability density is understood to be computed over all 
possible paths through the medium with an initial vertex distribution given by the nuclear overlap, 
i.e. the density of binary collisions. The detailed shape of P(AE,E) given a single path in the 
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medium can be computed in various models. In a recent comparison study [jUj], the medium 
modification of hadron spectra has been computed in three different frameworks with the averaging 
done over a 3-d hydrodynamical evolution model of the bulk medium. In each model, there is 
one free parameter which links the strength of the medium modification to the thermodynamical 
variables in the hydrodynamical model. This parameter is adjusted to data. The result in terms of 
the nuclear suppression factor Raa is shown in Fig. [l| 
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Figure 1: Left panel: Nuclear suppression factor Raa averaged over the angle with the reaction plane as 
a function of hadron momentum Pt in 200 AGeV Au-Au collisions in three different framworks compared 
with PHENIX data for two different centralities. Right panel: Raa calculated as a function of Pt in noncenttal 
200 AGeV Au-Au collisions assuming an impact parameter b = 7.5 fm for in plane (0 = 0) and out of plane 
(<j) = 90) emission where <p is the angle of the outgoing parton with the reaction plane (figures from JlT|]). 



It is clear from these results that all models describe both the Pj and the centrality dependence 
of the data well, although they exhibit some differences in the in plane vs. out of plane suppression 
which could be used to consttain the physics of the parton-medium interaction more. This similarity 
of the results, in spite of different physics assumptions entering the models, strongly emphasizes the 
need to go beyond the energy loss approximation. However, already leading hadron suppression 



allows some constraints on the physics of parton-medium interaction. In [20], using the energy 
loss approximation for both single hadron and back-to-back correlation suppression, it could be 
shown that the parametric pathlength dependence of the medium modification agrees with the idea 
that the dominant medium effect is induced gluon radiation, but does not agree with the idea that 
the parton interacts with the medium chiefly by elastic collisions. Tentatively, the contribution of 
elastic processes could be constrained to be about 10%. 



4. In-medium shower evolution and jets 

There are several advantages to studying jets instead of single inclusive hadron spectra, in 
spite of the experimental difficulty of identifying jets above the event background in a heavy-ion 
collision. First, jets reflect unbiased hard events in which the momentum flow is not constrained 
to be through a single parton, hence given the same shower momentum, they can be measured 
with higher statistics, essentially leading to a larger experimental reach in Pt- Furthermore, jet 
observables allow to observe the energy flow from leading to subleading partons and hence offer a 
more detailed window on medium effects than energy loss from a leading parton. Finally, shower 
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evolution is usually computed in momentum space, but the medium provides a 'meter stick' in 
terms of its finite extension and a 'clock' in terms of its finite lifetime to observe the evolution also 
in position space. 

The last property actually makes it necessary for any model of in-medium shower evolution 
to link the momentum space evolution outlined above with the propagation of partons in position 
space where the bulk medium is evolved. In the following, I outline the physics assumptions 



underlying the MC code YaJEM. A detailed description of the model can be found in \ \14[ |lq ]. In 
YaJEM the link to the spacetime dynamics is made by modelling the average time for a parton b to 
branch from parent a given the parton energies and virtualities based on the uncertainty relation as 

whereas the actual time in given branching is generated from the exponential branching prob- 
ability distribution 



P(% b ) = exp 



(4.2) 



For simplicity, all partons of a shower are propagated with this time information along an 
eikonal trajectory determined by the shower initiator. This amounts to neglecting the spread in 
transverse space when probing the medium. 

Currently, YaJEM models three different scenarios for the parton-medium interaction, two 
of which modify the kinematics of the propagating parton whereas the last modifies the branching 
probabilities at each vertex (thus, in the first two scenarios the energy of the shower is not the energy 
of the shower initiating parton as there is explicit energy transfer between shower and medium, 
whereas in the last scenario the energy in the shower is conserved). In the RAD scenario, the 
medium is assumed to cause an increase AQl in the virtuality of a parton a based on a local transport 
coefficient q(Q as given by the line integral 

AQl= f W d^q(Q. (4.3) 
This modification causes medium-induced radiation. Another modification, referred to as 



DRAG is a drag coefficient D(Q which is motivated by the result of AdS/CFT computations Q2T| ] 
of energy loss in a strongly coupled medium. This induces energy loss AE a for every propagating 
parton a according to the line integral 

AE a = / " "dCDpiQ. (4.4) 



The last scenario, called FMED, is included for comparison with the modelling done in [13 
. Here, the singular part of the splitting kernel in the evolution equations is enhanced by a factor 
(1 + fmed) where the assumption f mec i ~ J d£p(Q provides the link with the spacetime evolution 
in YaJEM. 

Fig. ^ left panel shows the medium-modified fragmentation function for a 20 GeV <i-quark for 
a path through the medium in all three scenarios characterized by the specific medium parameters, 
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Figure 2: Left panel: Medium-modified fragmentation function of a 20 GeV c/-quark computed in YaJEM 
for the vacuum and three different scenarios of parton-medium interaction (see text), adjusted such that the 
fragmentation functions agree around z ~ 0.7. Right panel: Nuclear modification factor Raa for central 200 
AGeV Au-Au collisions computed by averaging over a 3-d hydrodynamical model of the medium for three 
scenarios of parton-medium interaction. 



where these parameters have been selected for the purpose of comparison such that the fragmenta- 
tion functions agree around z ~ 0.7. It is clearly evident that the medium causes a depletion of the 
fragmentation function at high z. This is, in the language of full in-medium shower evolution, the 
equivalent of energy loss. In the right panel, Raa is computed by averaging over all possible paths 
through the medium. The medium-induced suppression is clearly visible, although the decreasing 
trend with Pj is somewhat puzzling. However, the same trend is also observed in other calculations 
of full in-medium shower evolution. 

Based on these results, it would seem difficult to distinguish the different scenarios. However, 
one can also look at other jet observables, such as the longitudinal momentum distribution inside 
the jet in terms of ^ = ln[l /z] which magnifies the low z region of the fragmentation function or 
the angular distribution of hadrons in the jet. This is shown in Fig. |3[ 




Figure 3: Left panel: dN/dE, for vacuum and three different scenarios of in-medium shower evolution (see 
text) Right panel: The angular distribution of hadrons in the shower with respect to the jet axis above 1 GeV 
momentum the vacuum and three different scenarios of in-medium shower evolution (see text). 
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For the two scenarios generating medium-induced radiation, RAD and FMED, a clear en- 
hancement in the high E, region over the vacuum result is visible in dN/d^, whereas this is not 
seen for the DRAG result. This corresponds to the presence or absence of multiple low Pj hadron 
production in the hadronization of the medium induced radiation. A measurement of this observ- 
able could distinguish between a radiative energy loss mechanism and a drag force which are able 
to produce a near-identical Raa- 

A similar case can be made with the angular distribution of hadrons, depicted in Fig. ||, right 
panel. Here, the radiative scenarios RAD and FMED lead to a widening of the jet whereas the 
DRAG scenario does not. If this could be measured, it could be a handle on the physics process 
underlying energy the parton-medium interaction. 

There is, however, an important caveat: Any experimental measurement of jets in a heavy- 
ion environment will have to make use of cuts to identify the jet. These cuts in turn will cause 
a systematic bias on the jet sample that is studied. In [15] it was shown that if the set of cuts is 
chosen too tightly, the resulting bias makes the medium modification of an observable (in the paper 
dN/dE, was considered) insignificant. 



5. Summary 

In this paper, I have presented arguments that jet observables contain information beyond what 
can be gained from a study of single inclusive hadron spectra and dihadron correlations in heavy- 
ion collisions. From arguments based on the uncertainty relation, one can deduce that the medium 
modification of a hard process chiefly affects the partonic evolution of a shower. For observables 
sensitive to the leading hadron of a jet only, the dominant kinematics is such that a single parton 
carries most of the momentum in the shower. As a consequence, the in-medium evolution of the 
shower can be modelled as energy loss from the leading parton. While computations based on this 
approximation are very successful when confronted with data, it remains difficult to distinguish dif- 
ferent physics assumptions underlying models. An exception to this is the pathlength dependence 
of energy loss which can be constrained well from the suppression of back-to-back correlations and 
which agrees with medium-induced radiation as the main mechanism of energy loss, but not with 
elastic collisions. 

Jet observables have the potential to reveal more about the nature of parton-medium interac- 
tion, as they trace the redistribution of energy in the medium. Since explicit energy-momentum 
conservation may be the most important factor in understanding this redistribution, the theoretical 
investigation of these observables is best made by MC simulations of the in-medium parton shower 
evolution. In this paper, I have outlined the physics underlying the MC shower code YaJEM. First 
results applying the code to single hadron observables or jet propertis find indeed the expected re- 
sult — while the nuclear suppression factor Raa is unable to distinguish between medium-induced 
raditation and a drag force as the effect causing the medium modification, jet observables such as 
the angular distribution of hadrons with respect to the jet axis or the longitudinal distribution of 
momenta in the jet show clear differences. 

There is now good reason to assume that jet observables will reveal important information 
about the microscopical properties of the medium, and Monte Carlo simulations of in-medium 
showers such as YaJEM will most likely be the appropriate tools to extract this information. 



8 



Medium-modified Jets 



Thorsten Renk 



References 

[1] M. Gyulassy and X. N. Wang, Nucl. Phys. B 420, (1994) 583; R. Baier, Y. L. Dokshitzer, 

A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 484, (1997) 265; B. G. Zakharov, JETP Lett. 
65, (1997) 615; U. A. Wiedemann, Nucl. Phys. B 588, (2000) 303; M. Gyulassy, P. Levai and I. Vitev, 
Nucl. Phys. B 594, (2001) 371; X. N. Wang and X. F. Guo, Nucl. Phys. A 696, (2001) 788. 

[2] M. Shimomura [PHENIX Collaboration], Nucl. Phys. A 774 (2006) 457. 

[3] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 76 (2007) 034904. 

[4] D. Magestro [STAR Collaboration], Nucl. Phys. A 774 (2006) 573; J. Adams et al. [STAR 
Collaboration], Phys. Rev. Lett. 97 (2006) 162301. 

[5] T. Renk, J. Ruppert, C. Nonaka and S. A. Bass, Phys. Rev. C 75 (2007) 031902. 

[6] T. Renk, Phys. Rev. C 74 (2006) 024903. 

[7] T. Renk and K. Eskola, Phys. Rev. C 75 (2007) 054910. 

[8] A. Majumder, C. Nonaka and S. A. Bass, Phys. Rev. C 76 (2007) 041902. 

[9] G. Y. Qin, J. Ruppert, S. Turbide, C. Gale, C. Nonaka and S. A. Bass, Phys. Rev. C 76 (2007) 064907. 
[10] H. Zhang, J. F. Owens, E. Wang and X. N. Wang, Phys. Rev. Lett. 98 (2007) 212301. 
[11] S. A. Bass, C. Gale, A. Majumder, C. Nonaka, G. Y. Qin, T. Renk and J. Ruppert, 0808.0908 [nucl-th]. 
[12] J. Putschke [STAR Collaboration], 0809.1419 [nucl-ex]. 

[13] K. Zapp, G. Ingelman, J. Rathsman, J. Stachel and U. A. Wiedemann, 0804.3568 [hep-ph]. 
[14] T. Renk, Phys. Rev. C 78 (2008) 034908. 
[15] T. Renk, 0808.1803 [hep-ph]. 
[16] T. Renk, 0901.2818 [hep-ph]. 

[17] N. Armesto, L. Cunqueiro and C. A. Salgado, 0809.4433 [hep-ph]. 

[18] M. Bengtsson and T. Sjostrand, Phys. Lett. B 185 (1987) 435; Nucl. Phys. B 289 (1987) 810; 
E. Norrbin and T. Sjostrand, Nucl. Phys. B 603 (2001) 297. 

[19] T. Sjostrand, Comput. Phys. Commun. 82 (1994) 74. 

[20] T. Renk, Phys. Rev. C 76 (2007) 064905. 

[21] S. S. Gubser, Phys. Rev. D 74 (2006) 126005; C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and 
L. G. Yaffe, JHEP 0607 (2006) 013. 

[22] N. Borghini and U. A. Wiedemann, hep-ph/0506218. 



9 



